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STATEMENT  OF  WORK 


The  purpose  of  this  research  project  is  to  study  the  time  domain  response  of  electro¬ 
magnetic  wave  radiation,  transmission,  and  coupling  in  multilayered  media.  The  following 
problems  are  pursued:  (1)  extensions  and  modifications  to  the  double  deformation  tech¬ 
nique;  (2)  propagation  in  nonconventional  transmission  structures;  (3)  signal  distortion  at 
discontinuities;  (4)  the  effects  of  anisotropic  material  and  nonlinear  loads;  (5)  limitation  of 
quasi-TEM  approximation.  We  shall  emphasize  and  seek  to  refine  a  powerful  transform- 
domain  formulation,  the  double  deformation  technique  in  order  to  have  a  unified  way  of 
interpreting  the  results.  Yet  other  techniques  such  as  the  space-time  domain  integral  equa¬ 
tion  method,  the  transmission  line  matrix  method,  the  method  of  characteristics,  and  the 
method  of  moments  axe  also  to  be  applied  to  different  problems  as  demanded  by  efficiency 
or  ease  of  formulation.  The  research  results  can  be  applied  to  computer-aided  design  of 
high-speed  microelectronic  integrated  circuits,  as  well  as  to  time-domain  geophysics  sub¬ 
surface  probing,  and  active  remote  sensing  with  transient  radar  pulses. 


SUMMARY  OP  RESEARCH  FINDINGS 


The  guidance  and  leakage  properties  of  single  and  coupled  dielectric  strip  waveg¬ 
uides  are  analyzed  using  the  dyadic  Green’s  function  and  integral  equation  formulation. 
Galerkin’s  method  is  used  to  solve  the  integral  equation  for  the  dispersion  relation.  The 
effects  of  the  geometrical  and  the  electrical  parameters  on  the  dispersion  relation  are  inves¬ 
tigated.  A  method  to  predict  the  occurrence  of  leakage  is  proposed.  The  properties  of  the 
even  and  the  odd  leaky  modes  are  also  investigated.  Results  are  compared  with  previous 
analysis  and  shown  to  be  in  good  agreement. 

The  input  impedance  of  a  microstrip  antenna  consisting  of  two  circular  microstrip 
disks  in  a  stacked  configuration  driven  by  a  coaxial  probe  is  investigated.  A  rigorous  anal¬ 
ysis  is  performed  using  a  dyadic  Green’s  function  formulation  where  the  mixed  boundary 
value  problem  is  reduced  to  a  set  of  coupled  vector  integral  equations  using  the  vector 
Hankel  transform.  Galerkin’s  method  is  employed  in  the  spectral  domain  where  two  sets 
of  disk  current  expansions  are  used.  One  set  is  based  on  the  complete  set  of  orthogonal 
modes  of  the  magnetic  cavity,  and  the  other  employs  Chebyshev  polynomials  with  the 
proper  edge  condition  for  the  disk  currents.  An  additional  term  is  added  to  the  disk  cur¬ 
rent  expansion  to  properly  model  the  current  in  the  vicinity  of  the  probe/disk  junction. 
The  input  impedance  of  the  stacked  microstrip  antenna  including  the  probe  self-impedance 
is  calculated  as  a  function  of  the  layered  parameters  and  the  ratio  of  the  two  disk  radii. 
Disk  current  distributions  and  radiation  patterns  are  also  presented.  The  calculate  results 
are  compared  with  experimental  data  and  shown  to  be  in  good  agreement. 

The  frequency-dependent  resistance  and  inductance  of  uniform  transmission  lines 
are  calculated  with  a  hybrid  technique,  which  is  a  combination  of  a  cross-section  coupled 
circuit  method  and  a  surface  integral  equation  approach.  The  coupled  circuit  approach 
is  most  applicable  for  low-frequency  calculations,  while  the  integral  equation  approach  is 
best  for  high  frequencies.  The  low-frequency  method  consists  of  subdividing  the  cross- 
section  of  each  conductor  into  triangular  filaments,  each  with  an  assumed  uniform  current 
distribution.  The  resistance  and  mutual  inductance  between  the  filaments  are  calculated, 
and  a  matrix  is  inverted  to  give  the  overall  resistance  and  inductance  of  the  conductors. 
The  high-frequency  method  expresses  the  resistance  and  inductance  of  each  conductor  in 
terms  of  the  current  at  the  surface  of  that  conductor  and  the  derivative  of  that  current 
normal  to  the  surface.  A  coupled  integral  equation  is  then  derived  to  relate  these  quantities 
through  the  diffusion  equation  inside  the  conductors  and  Laplace’s  equation  outside.  The 
method  of  moments  with  pulse  basis  functions  is  used  to  solve  the  integral  equations.  An 


interpolation  between  the  results  of  these  two  methods  gives  very  good  results  over  the 
entire  frequency  range,  even  when  few  basis  functions  are  used. 

Because  the  effects  of  diffraction  during  proximity-print  x-ray  lithography  are  of 
critical  importance,  a  number  of  previous  researchers  have  attempted  to  calculate  the 
diffraction  patterns  and  minimum  achievable  feature  sizes  as  a  function  of  wavelength  and 
gap.  Work  to  date  has  assumed  that  scalar  diffraction  theory  is  applicable-as  calculated,  for 
example,  by  the  Rayleigh- Sommerfeld  formulation-and  that  Kirchhoff  boundary  conditions 
can  be  applied.  Kirchhoff  boundary  conditions  assume  that  the  fields  (amplitude  and 
phase)  are  constant  in  the  open  regions  between  absorbers,  and  a  different  constant  in 
regions  just  under  the  absorbers  (i.e.,  that  there  are  no  fringing  fields).  An  x-ray  absorber 
is,  however,  best  described  as  a  lossy  dielectric  that  is  tens  or  hundreds  of  wavelengths 
tall,  and  hence  Kirchhoff  boundary  conditions  are  unsuitable.  We  found  out  that  the  use 
of  Kirchhoff  boundary  conditions  introduces  unphysically  high  spatial  frequencies  into  the 
diffracted  fields.  The  suppression  of  these  frequencies-which  occurs  naturally  without  the 
need  to  introduce  an  extended  source  or  broad  spectrum-improves  exposure  latitude  for 
mask  features  near  0.1  fim  and  below. 

The  electromagnetic  radiation  from  a  VLSI  chip  package  and  heatsink  structure 
is  analysed  by  means  of  the  finite-difference  time-domain  (FD-TD)  method.  The  FD- 
TD  algorithm  implemented  incorporates  a  multi-zone  gridding  scheme  to  accommodate 
fine  grid  cells  in  the  vicinity  of  the  heatsink  and  package  cavity  and  sparse  gridding  in 
the  remainder  of  the  computational  domain.  The  issues  pertaining  to  the  effects  of  the 
heatsink  in  influencing  the  overall  radiating  capacity  of  the  configuration  are  addressed. 
Analyses  are  facilitated  by  using  simplified  heatsink  models  and  by  using  dipole  elements  as 
sources  of  electromagnetic  energy  to  model  the  VLSI  chip.  The  potential  for  enhancement 
of  spurious  emissions  by  the  heatsink  structure  is  illustrated.  For  heatsinks  of  typical 
dimensions,  resonance  is  possible  within  the  low  gigahertz  frequency  range.  The  potential 
exploitation  of  the  heatsink  as  an  emissions  shield  by  appropriate  implementation  schemes 
is  discussed  and  evaluated. 

A  method  for  the  calculation  of  the  current  distribution,  resistance,  and  inductance 
matrices  for  a  system  of  coupled  superconducting  transmission  lines  having  finite  rectan¬ 
gular  cross  section  is  derived.  These  calculations  allow  accurate  characterization  of  both 
high-Te  and  low-Te  superconducting  strip  transmission  lines.  For  a  single  stripline  geome¬ 
try  with  finite  ground  planes,  the  current  distribution,  resistance,  inductance,  and  kinetic 
inductance  are  calculated  sis  a  function  of  the  penetration  depth  for  various  film  thickness. 
These  calculations  are  then  used  to  determine  the  penetration  depth  for  Nb,NbN,  and 
YBajCtttOr-,  superconducting  thin  films  from  the  measured  temperature  dependence  of 
the  resonant  frequency  of  a  stripline  resonator.  The  calculations  are  also  used  to  convert 
measured  temperature  dependence  of  the  quality  factor  to  the  intrinsic  surface  resistance 
as  a  function  of  temperature  for  a  Nb  stripline  resonator. 


A  general  spectral  domain  formulation  to  the  problem  of  radiation  of  arbitrary 
distribution  of  sources  embedded  in  a  horizontally  stratified  arbitrary  magnetized  linear 
plasma  is  presented.  The  fields  are  obtained  in  terms  of  electric  and  magnetic  type  dyadic 
Green’s  functions.  The  formulation  is  considerably  simplified  by  using  the  kDB  system 
of  coordinates  in  conjunction  with  the  Fourier  transform.  The  distributional  singular 
behavior  of  the  various  dyadic  Green’s  functions  in  the  source  region  is  investigated  and 
taken  into  account  by  extracting  the  delta  function  singularities.  Finally,  the  fields  in  any 
arbitrary  layer  are  obtained  in  terms  of  appropriately  defined  global  upward  and  downward 
reflection  and  transmission  matrices. 

A  spectral  domain  dyadic  Green’s  function  formulation  defining  the  fields  inside 
a  multilayer  chiral  medium  due  to  arbitrary  distribution  of  sources  is  presented.  The 
constitutive  parameters  and  the  chirality  of  each  layer  are  assumed  to  be  different.  The 
fields  are  obtained  in  terms  of  electric  and  magnetic  type  dyadic  Green’s  functions.  The 
singular  behavior  of  these  dyadic  Green’s  functions  in  the  source  region  is  taken  into 
account  by  extracting  the  delta  function  singularities.  The  fields  in  any  layer  are  obtained 
in  terms  of  appropriately  defined  global  reflection  and  transmission  matrices. 

The  frequency-dependent  resistance  and  inductance  of  uniform  transmission  lines 
are  calculated  using  a  hybrid  technique,  which  is  a  combination  of  a  cross-section  finite 
element  method  and  a  surface  integral  equation  approach.  The  finite  element  approach  is 
most  applicable  for  low-frequency  calculations,  while  the  integral  equation  approach  is  best 
for  high  frequencies.  An  interpolation  between  the  results  of  these  two  methods  gives  very 
good  results  over  the  entire  frequency  range,  even  when  few  basis  functions  are  used.  Using 
this  method,  a  potential  CAD  tool  for  the  calculation  of  transmission  line  parameters  is 
developed. 

A  direct  three-dimensional  finite- difference  time-domain  (FDTD)  method  is  applied 
to  the  full-wave  analysis  of  various  microstrip  structures.  The  method  is  shown  to  be  an 
efficient  tool  for  modeling  complicated  microstrip  circuit  components  as  well  as  microstrip 
antennas.  From  the  time-domain  results,  the  input  impedance  of  a  line-fed  rectangular 
patch  antenna  and  the  frequency-dependent  scattering  parameters  of  a  low-pass  filter  and  a 
branch  line  coupler  are  calculated.  These  circuits  are  fabricated  and  the  measurements  are 
compared  with  the  FDTD  results  and  shown  to  be  in  good  agreement.  A  general  purpose 
time-domain  numerical  algorithm  for  modelling  three-dimensional  microstrip  structures  is 
developed. 

The  pseudo-differential  operator  approach  is  employed  to  derive  absorbing  bound¬ 
ary  conditions  for  both  circular  and  elliptical  outer  boundaries.  The  pseudo-differential 
operator  approach  employed  by  Engquist  and  Majda  is  modified  to  derive  improved  ab¬ 
sorbing  boundary  conditions.  In  the  case  of  circular  outer  boundaries,  the  modified  pseudo¬ 
differential  operator  approach  leads  to  a  condition  equivalent  to  that  of  Bayliss  and  Turkel’s 


second-order  condition.  The  modified  pseudo-differential  operator  is  then  used  to  derive 
the  second-order  absorbing  boundary  condition  for  elliptical  outer  boundaries.  The  effec¬ 
tiveness  of  the  second-order  absorbing  boundary  condition  on  elliptical  outer  boundary 
is  illustrated  by  calculating  scattered  fields  from  various  objects.  It  is  shown  that  for 
elongated  scatterers,  the  elliptical  outer  boundary  can  be  used  to  reduce  the  s:ze  of  the 
computational  domain. 

Quasi-TEM  approximation  is  applied  to  the  analysis  of  coupled  lossy  microstrip 
lines  with  finite  thickness  embedded  in  a  horizontally  stratified  medium.  A  scalar  Green’s 
function  in  the  spectral  domain  is  used  to  obtain  a  set  of  coupled  integral  equations  for 
the  surface  charge  distribution.  The  method  of  moments  is  then  used  to  find  the  charge 
distribution  and  hence  the  capacitance  matrix  of  the  microstrip  lines.  The  inductance  and 
the  conductance  matrices  are  obtained  by  using  the  duality  between  the  magnetostatic 
problem,  the  current  field  problem,  and  the  electrostatic  problem.  The  resistance  matrix 
is  obtained  by  a  perturbation  method.  A  multiconductcr  transmission  line  analysis  is 
derived  by  using  the  capacitance,  the  inductance,  the  conductance,  and  the  resistance 
matrices.  The  transient  response  is  obtained  by  using  the  Fourier  transform. 

An  inversion  algorithm  based  on  a  recently  developed  inversion  method  referred  to 
as  the  Renormalized  Source-Type  Integral  Equation  approach  is  developed.  The  objective 
of  this  method  is  to  overcome  some  of  the  limitations  and  difficulties  of  the  iterative  Born 
technique.  It  recasts  the  inversion,  which  is  nonlinear  in  nature,  in  terms  of  the  solution 
of  a  set  of  linear  equations;  however,  the  final  inversion  equation  is  still  nonlinear.  The 
derived  inversion  equation  is  an  exact  equation  which  sums  up  the  iterative  Neuman  (or 
Born)  series  in  a  closed  form  and;  thus,  is  a  valid  representation  even  in  the  case  when  the 
Born  series  diverges. 

The  coupled-wave  theory  is  generalized  to  analyze  the  diffraction  of  waves  by  chiral 
gratings  for  arbitrary  angle  of  incidence  and  polarizations.  Numerical  results  are  obtained 
for  the  Stokes  parameters  of  diffracted  Floquet  modes  versus  the  thickness  of  chiral  grat¬ 
ings  with  various  chiralities.  Both  horizontal  and  vertical  incidences  are  considered.  The 
diffracted  waves  from  chiral  gratings  are  in  general  elliptically  polarized;  and  at  some  par¬ 
ticular  instances,  it  is  possible  for  chiral  gratings  to  convert  a  linearly  polarized  incident 
field  into  two  nearly  circularly  polarized  Floquet  modes  propagating  in  different  directions. 

A  theoretical  analysis  of  a  circular  microstrip  antenna  consisting  of  two  disks  in  a 
stacked  configuration  driven  by  a  coaxial  probe  is  presented.  The  resonant  frequencies, 
input  impedance,  radiation  fields,  and  receiving  and  scattering  characteristics  are  inves¬ 
tigated.  For  each  of  the  problems  presented,  a  rigorous  analysis  is  performed  using  a 
dyadic  Green’s  function  formulation.  Employing  the  vector  Hankel  transform,  the  prob¬ 
lem  is  reduced  to  a  set  of  coupled  vector  integral  equations  and  solved  using  Galerkin’s 
method  in  the  spectral  domain.  The  scattering,  receiving,  and  transmitting  characteristics 


of  an  infinite  array  of  probe-fed  stacked  circular  microstrip  antennas  are  also  investigated. 
The  input  impedance,  scanning  performance,  and  cross  polarization  levels  are  thoroughly 
investigated. 
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Input  Impedance  of  a  Probe-Fed  Stacked 
Circular  Microstrip  Antenna 
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Absrroct—Tbt  iapat  hnpwhace  of  a  mJcrostrip  antcaaa  coexisting  of 
two  drcatar  microstrip  disks  in  a  stacked  configuration  drive*  by  a 
coaxial  probe  is  iarestigaied.  A  rigorous  analysis  is  performed  a«iag  a 
dyadic  Green’s  function  fomMiatioa  where  the  mixed  boundary  mine 
problem  is  reduced  to  a  set  of  coopted  vector  integral  equations  asing 
the  vector  Hank  el  transform.  Galerkin's  method  is  employed  h  the 
spectral  domain  where  two  sets  of  disk  cwnreat  expansions  art  used.  One 
set  is  based  on  the  complete  set  of  orthogonal  modes  of  the  magnetic 
cavity,  aad  the  other  emptoyt  Chebyshev  polynomials  with  the  proper 
edge  condition  for  the  disk  currents.  An  additional  term  is  added  to  the 
disk  earnest  expansion  to  properly  model  the  current  is  the  vicinity  of 
the  probe/disk  Junction.  The  input  impedance  of  the  stacked  microstrip 
antenna  including  the  probe  self-impedance  is  calculated  as  a  function 
of  the  layered  parameters  aad  the  ratio  of  the  two  disk  radii.  Disk 
current  distributions  and  radiation  patterns  are  also  presented.  The 
calculated  results  are  compared  with  experimental  data  and  shown  to  be 
hs  good  agreement 


I.  Introduction 

CONVENTIONAL  microstrip  antennas,  consisting  of  a  sin¬ 
gle  conducting  patch  on  a  grounded  dielectric  substrate, 
have  received  much  attention  in  recent  years  [1]  due  to  their 
many  advantages,  including  low  profile,  light  weight,  and  easy 
integration  with  printed  circuits.  However,  due  to  their  resonant 
behavior,  they  radiate  efficiently  only  over  a  narrow  band  of 
frequencies,  with  bandwidths  typically  only  a  few  percent  [1J. 
While  maintaining  the  advantages  of  conventional  single  patch 
microstrip  antennas,  microstrip  antennas  of  stacked  configura¬ 
tions,  consisting  of  one  or  more  conducting  patches  parasitically 
coupled  to  a  driven  patch,  overcome  the  inherent  narrow  band¬ 
width  limitation  by  introducing  additional  resonances  in  the 
frequency  range  of  operation,  achieving  bandwidths  up  to 
10-20%.  In  addition,  stacked  microstrip  configurations  have 
achieved  higher  gains  and  offer  the  possibility  of  dual  frequency 
operation. 

Experimental  work  with  multilayered  microstrip  elements  has 
been  abundant  [2] -[9].  However,  to  date,  theoretical  work  has 
been  relatively  limited,  where  the  study  of  resonant  frequencies, 
modes  and  radiation  patterns  have  been  investigated  [10]-(13]. 
Recently,  the  finite-difference  time-domain  technique  was  ap¬ 
plied  to  stacked  rectangular  microstrip  patch  configurations  [14]. 
There  is  little  or  no  theoretical  analysis  of  the  input  impedance 
of  coaxial  probe-fed  stacked  circular  microstrip  patches.  How¬ 
ever,  the  input  impedance  for  conventional  single-layer  coaxial 
probe-fed  microstrip  antennas  of  circular,  rectangular,  annular 
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Fig.  1 .  Stacked  microstrip  antenna  configurations. 


ring,  and  elliptic  geometries  has  been  investigated  by  many 
authors  [15]-[19],  The  impedance  parameters  of  two  planar 
coupled  microstrip  patches  have  also  been  studied  [19],  [20]. 

In  the  calculation  of  the  input  impedance  of  probe-driven 
microstrip  antennas  on  thin  substrates,  the  effect  of  the  probe 
results  in  an  additional  inductive  component  to  the  input 
impedance.  This  probe  inductance  has  been  accounted  for  by 
several  authors  through  use  of  a  simple  formula  [19],  [21].  In 
more  rigorous  methods  to  include  the  effects  of  the  probe,  an 
“attachment  mode"  in  the  disk  current  expansion  is  used  to 
account  for  the  singular  behavior  of  the  disk  current  in  the 
vicinity  of  the  probe,  ensure  continuity  of  the  current  at  the 
probe/disk  junction,  and  speed  up  the  convergence  of  the  solu¬ 
tion.  An  "attachment  mode”  which  represented  the  disk  current 
of  a  lossy  magnetic  cavity  driven  by  a  uniform  cylindrical  probe 
current  was  introduced  in  [16].  More  recently,  similar  and  other 
“attachment  modes,”  with  the  1  /p  dependence  in  the  vicinity  of 
the  probe  and  the  appropriate  boundary  condition  on  normal 
current,  defined  over  the  entire  disk  or  locally  over  a  portion  of 
the  disk,  have  also  been  used  [22]— [25].  In  a  different  approach, 
the  effects  of  the  probe  were  accounted  for  by  expanding  the 
currents  on  the  disk  and  probe  in  terms  of  the  modes  of  a 
cylindrical  magnetic  cavity  satisfying  boundary  conditions  on  the 
eccentrically  located  probe  [26]. 

Considered  here  is  a  microstrip  antenna  consisting  of  two 
circular  microstrip  disks  in  a  stacked  configuration  driven  by  a 
coaxial  probe.  The  two  stacked  configurations  shown  in  Figs. 
1(a)  and  1(b),  denoted  configurations  A  and  B,  respectively,  are 
investigated.  The  disks  are  assumed  to  be  infinitesimally  thin  and 
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perfectly  conducting  and  the  substrates  are  taken  to  be  infinite  in 
extent.  A  rigorous  analysis  of  the  two  stacked  circular  disks  in  a 
layered  medium  is  performed  using  a  dyadic  Green’s  function 
formulation.  Using  the  vector  Hankel  transform,  the  mixed 
boundary  value  problem  is  reduced  to  a  set  of  coupled  vector 
integral  equations  and  solved  by  employing  Galerldn's  method  in 
the  spectral  domain.  Two  solutions  using  two  different  basis  sets 
to  expand  the  unknown  disk  currents  are  developed.  The  first  set 
of  basis  functions  used  are  the  complete  set  of  transverse  mag¬ 
netic  (TM)  and  transverse  electric  (TE)  modes  of  a  cylindrical 
cavity  with  magnetic  side  walls.  The  second  set  of  basis  func¬ 
tions  used  employ  Chebyshev  polynomials  and  enforce  the  cur¬ 
rent  edge  condition.  An  additional  term  in  the  current  expansion 
is  taken  to  account  for  the  singular  nature  of  the  current  on  the 
disk  in  the  vicinity  of  the  probe  and  to  ensure  continuity  of 
current  at  the  junction.  This  term,  the  “attachment  mode,"  is 
taken  to  be  the  disk  current  of  a  magnetic  cavity  under  a  uniform 
cylindrical  current  excitation.  It  is  shown  here  explicitly  that 
continuity  of  the  current  at  the  probe /disk  junction  must  be 
enforced  to  rigorously  include  the  probe  self-impedance.  The 
convergence  of  the  results  is  investigated  and  ensured  by  using  a 
proper  number  of  basis  functions.  The  input  impedance  of  the 
stacked  microstrip  antenna  is  calculated  for  different  configura¬ 
tions  of  substrate  parameters  and  disk  radii.  Disk  current  distri¬ 
butions  and  radiation  patterns  are  also  presented.  Finally,  the 
results  are  compared  with  experimental  data  and  shown  to  be  in 
good  agreement.  Throughout  the  analysis,  the  exp(-iwt)  time 
dependence  is  used  and  suppressed. 

II.  Dyadic  Green’s  Function  and  Integral 
Equation  Formulation 

For  a  general  formulatiori  which  applies  to  both  configurations 
A  and  B  of  Fig.  1,  we  consider  two  coaxial,  circular  perfectly 
conducting  disks,  of  radii  a,  and  a2,  carrying  current  distribu¬ 
tions  Jj(r)  =  Kj(p)h(z  -  Zj)  where  j  =  (1,2)  and  $(•)  is  the 
Dirac  delta  function.  Configuration  A  is  obtained  when  z\  -  </, 
and  z'2  -  d2  and  configuration  B  results  when  zj  *  d0  and 
t2  »  d2. 

Using  the  induced  EMF  method  [27],  a  stationary  formula  for 
the  input  impedance  is  obtained  as 


equations  for  the  disk  currents 
[£(p,  z  =  z;)]r=  £  eim*  ja‘ dkfi  kpJm(kpfi) 
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where  j  =  1,2,  and  kp  is  the  transverse  wavenumber  satisfying 
the  dispersion  relation 

+  ku  -  */  =  “V,  (5) 

in  each  region  /.  In  (3),  z  and  z'  correspond  to  the  longitudinal 
positions  of  observer  and  source,  respectively.  k^>( kp)  and 
•^m\kp)  are  the  vector  Hankel  transforms  of  the  two  disk 
currents  K£\p)  and  K(*\p),  respectively,  defined  by 

=  r dp pJl{k,p)  ■  E$(p)  (6) 

Jo 

where  is  the  Fourier  coefficient 

K}?(p)-T-  [2'd*e-‘”»Kj(p)  (7) 
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and  Jm(kpp)  is  the  kernel  of  the  vector  Hankel  transform 
(VHT)  [28]  given  by 


z“  =  "  F  ///  dVE{f)  ‘  (1) 

where  is  die  current  distribution  on  the  probe  and  £  is  the 
total  electric  field  due  to  the  probe  current  and  induced  disk 
currents. 

The  current  on  the  probe,  of  radius  R  and  at  the  position 
p0  m  (p0,  is  taken  to  be  uniform  and  is  given  by 

d2<  z  <  d2  (2) 

with  local  coordinates  defined  as  pp  *  p  -  p0  ■  (pp,  <bp)- 
Using  a  dyadic  Green’s  function  formulation  in  cylindrical 
coordinates  for  horizontally  stratified  media  [11],  we  obtain 
expressions  for  the  transverse  components  of  the  electric  fields 
due  to  the  disk  and  probe  current  distributions.  Boundary  condi¬ 
tions  require  that  the  transverse  components  of  the  electric  field 
vanish  on  the  perfectly  conducting  disks  and  the  currents  vanish 
off  the  disks,  to  give  the  following  set  of  coupled  integral 


^m(kPP) 


J'm(k„p)  T7JAk,p) 
k„P 

im 

77-UM  Jm{k„p) 

k„p 


(8) 


Jm(  •)  is  the  Bessel  function  of  the  first  kind  of  order  m  and  the 
prime  denotes  differentiation  with  respect  to  the  argument. 

Jm(kPP)  is  the  complex  conjugate  transpose  of  Jm{kBp). 

In  the  last  term  of  (3),  pm(kfi)  is  associated  with  the  probe 
current  and  is  given  by 


Pm(kp) 

where 


Pm{k,) 

0 


(9) 


£«(*,)-  -T-TrJm(kl,Po)MkpX)*-im*0-  0°) 

LT  *3jt 

The  matrix  f^(kt,  z)  includes  the  effects  of  the  stratified 
medium  when  relating  the  probe  current  to  the  transverse  elec- 
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trie  fields  and  is  defined  as 

6s  (**.«) 


f™  0 
0  0 


(U) 


It  is  dear  that  the  assumed  probe  current  excites  TM  modes 
only.  The  matrices  z,  zO  with  l,j  =  (1,2)  in  (3)  in¬ 

clude  the  effects  of  the  stratified  medium  when  relating  the  disk 
currents  to  transverse  electric  fields  and  are  of  the  form 


Z,  Z  )  = 


£TM 

«/.> 


0 

t  TE 


(12) 


for  m  =  0, ±  1,±  2,---,  n  =  1,2,---,  and  p  =  1,2, 
p)  correspond  to  the  TM  cavity  modes  and  *^(p)  corre¬ 
spond  to  the  TE  cavity  modes.  The  constants  fimn  and  amp 
correspond  to  the  nth  and  pth  zeros  of  J'm(0mn)  =  0  and 
Jm( u„p)  =  0,  respectively.  The  VHT  of  these  basis  functions  is 


&£(*,)  =  ^,4(^„) 


-  A: 


rma 


r-U*„«>) 


The  expressions  for  {™(Arp,  z),  z,  zO.  and 

(™(Arp,  z,  zO  are  given  in  the  Appendix. 

HI.  Galerkin’s  Method 

Galeririn’s  method  is  employed  to  solve  the  coupled  vector 
integral  equations  of  (3)  and  (4).  The  currents  on  the  circular 
disks  are  expanded  in  terms  of  a  set  of  basis  functions 

*i!’(p)  =  Z&M  +  f«»p)  (13a) 

n  p 

*?(p)  =  £«2US(p)  +  I«(p)  +  *£!.„(*)• 

(13b) 

N  and  P  correspond  to  the  number  basis  functions  ^mj(p)  and 
<t>mJ(p),  respeciively,  taken  for  the  upper  disk_and  R  and  S 
correspond  to  those  taken  for  the  lower  disk.  K%'m(.p)  is  the 
“attachment  mode.” 

The  corresponding  VHT  of  the  currents  is  given  by 

sSW  =  +  !>»(*,)  04a) 

n  p 

*£>(*,)  =  tfirWAK)  +  it&PSAk,)  +  .»(*,)• 

(14b) 

i4.  7M  and  TE  Modes  of  Cylindrical  Cavities  with 
Magnetic  Side  Walls 

One  set  of  basis  functions  taken  are  those  currents  associated 
with  the  complete  orthogonal  set  of  TM  and  TE  modes  of  a 
cylindrical  cavity  of  radius  at  (j  =  1,2)  with  magnetic  side 
walls  and  electric  top  and  bottom  walls.  These  current  modes  are 
given  by 


<(*,)  - 


~  (amp/fly) 


0 


(16a) 


(16b) 


B.  Chebyshev  Polynomial  Expansion  with  Edge 
Condition 

The  second  set  of  basis  functions  taken  includes  the  edge 
condition  for  the  disk  currents  and  is  taken  to  be  [23] 


%2(p) 

3fi(p) 


|  pT„(p/ay)  y/l  -  p'/a)  ,  for  p  <  Oj 
|0,  for  p  >  Oj 


4T„(p/ Oj)/ \/\  -  p2/a*  ,  forpCfl, 
0,  for  p>  Oj 


(17a) 


(17b) 


for  m  *  0,±  1,±  2,  •  •  •  and  n  *  0, 1,2,  ••  •  T„{x)  is  the 
Chebyshev  polynomial  [29]  and  satisfies  the  recursion  formula 
rB+I(x)  -  2 xTn(x)  +  r„_,(x)  «  Q  with  T0(x)  =  1  and 

F,(  X)  =  x.  The  term  yj  1  -  p7  /  a*  provides  for  the  proper 
singular  edge  behavior  for  the  azimuthally  directed  current  and 
the  zero  edge  condition  for  the  normally  directed  current.  Since 
the  current  basis  functions  must  have  continuous  current  distri¬ 
butions  on  the  disk,  the  mode  index  m  and  the  Chebyshev 
polynomial  index  n  may  not  be  both  even  or  both  odd  when 
performing  the  current  expansion. 

The  VHT  of  the  above  basis  functions  is  given  by 


«£(*„) 


£■*(*,)  -  ^n\k„) 

-”*J&\k0) 
lUXk,)  ~  mjtiXk,) 


(18a) 

(18b) 
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/*«>.<  ±1 

— J-Jm{emnplaj ) 
PmnP 

7  (15a) 

mn  8  yi  2 

0, 

for  p  >  Oj 

’/;m  +  »i  +  2 

r 

T 

1  m  ■ f 

-(^1 


-ima 


LJm{^mpOlai) 


&mpP 

J'm(ampP/aj ) 
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,  for  p  <  fly 


for  p  >  fly 


/2)  Am-n-2)/l(yj  /2) 
'  |  A’n*n-7),r2(yj/2) 


(15b) 


'(m-n-f  2)f2 


( yj/2 ) 
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*  Jj  *• 

+  ^  J(  m + 1»)/2  (  JO  /2  )  -  «)/2  (  .P/  /2) 

~  ^(in+ji-O/lf  V2)  ,+2>/2(^/2)]  (19b) 

■*"'^(m+i»-2)/2(^’y/2)^(m-ii)/2(^’>/2)J  (19c) 
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In  the  above  expressions,  when  m  and  n  are  not  both  even  or 
both  odd,  the  Bessel  functions  are  of  half-integer  order 

I30J. 

C.  Attachment  Mode 

The  “attachment  mode"  term  in  the  current  expansion  is 
taken  to  approximate  the  rapidly  varying  currents  in  the  vicinity 
of  the  probe/disk  junction,  ensure  continuity  of  the  current,  and 
speed  up  the  convergence  of  the  solution.  This  term  is  taken  as 
the  disk  current  of  a  magnetic  cavity  of  radius  a2  due  to  a 
uniform  cylindrical  current  source  of  radius  R  positioned  at  p0 
and  given  by 

£  e"*— 0 

LT  m  *  -  ao  J  o  *3 z 


■J0(k,R)Jm{k,p  o) 


■£(M 

im 
kcP 


•UM 


ik2I 


+~r  £ 

**  m*  - 


■^m(^3a2) 


im 

*iP 


JJ*3  P) 


Ospst, 


(20) 


where  H^(-)  is  the  Hankei  function  of  the  first  land  of  order 
m.  The  first  term  in  (20)  is  the  current  induced  on  infinite 
parallel  conducting  planes  by  a  uniform  cylindrical  current.  The 
second  term  is  a  homogeneous  solutionjo  the  wave  equation 
added  to  satisfy  the  boundary  condition  //^(p  =  a2)  *  0,  pro¬ 
viding  for  vanished  normal  current  at  the  edge  of  the  disk. 

The  VHT  of  the  above  attachment  mode  current  distribution 
has  a  closed  form  analytic  expression  given  by 

*£-(*.)  - 

/0(*3*)'m(*3P0) 


2t 


^«(*3«2) 


kiz 

im 


*3«2*„ 


Uk.a2) 
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D.  Matrix  Equation 

Substituting  the  current  expansion  of  (14)  into  (3),  and  apply¬ 
ing  Parse val’s  theorem,  we  obtain  a  system  of  N  +  P  +  R  +  S 
linear  algebraic  equations  for  each  mode  m  which  may  be 
written  in  matrix  form 


where 


A  m  "  Cm  ~ 

3m 

NxN 

f 

l  NxP 

MUST 

Ml 

PxN 
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r  xi*(,V2> 
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and 


cm  ~ 
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Ml 
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(24) 


dm  = 


K] 

Nxl 

K”1 

Px  I 

K] 

Rx\ 

K 

Sxl  _ 

(25) 


Each  element  of  the  submatrices  of  A  m  is  given  by 


n mnp 


-  /V^(*,)  •  ■  MM 

J  0 

(26) 

where  ylJ^kp)  and  x^jkj  represent  either  #£(*„)  or 
$mn(k.).  Each  element  of  the  excitation  matrix  dm  is  given  by 


<nl,) 

umn 


-  /  ■  ijTik^z’,)  ■  Pm{ke) 

Jo 
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Jo 

(27) 


IV.  Input  Impedance 

Once  the  induced  current  distribution  on  the  microstrip  disks 
due  to  the  coaxial  probe  excitation  is  solved  for,  the  input 
impedance  of  the  stacked  microstrip  antenna  may  be  calculated. 
Applying  (1),  the  input  impedance  for  the  stacked  microstrip 
antenna  is  given  by 

-i  [d2dzJ2rd^rdpppp{[Ey,(f) ]t 

1  Jd3  Jo  Jo 

+  [««(')].+  [*-(')].»  •  27Ri{fi’~  R)  (28) 


(21) 
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where  E3J.r)  is  the  electric  field  due  to  disk  current  j  and 
E^f (r)  is  the  electric  field  due  to  the  probe  current.  After 
integration  over  the  cylindrical  probe  surface  and  some  manipu¬ 
lation,  we  arrive  at 

Z»=-7T  £  /  MPk,p*(k,){v™(kfi) 

i  m  «  -  ob  ^  o 

•(**!’(*,)],  +  nJ3(*,)[*J?(*,)U  + 

*-77  £  /  dkfikpP*(kt)L™(kp) 

■  p 

r 

+bB(*,)e«s[©m  J +  z<*n + 


waveguide.  In  the  small  k3R  limit,  this  term  reduces  to 


to  ■Tkth>Jo(k3R)Ho)(kiR) 

r\R-*0  4 


Ha  iir  ii « 

-  i—k0h3—  In  (*3R)  =  16O  k0h3—  In  {k3R) 

2*  Mo  Mo 

which  is  the  formula  used  by  some  authors  as  the  probe  reac¬ 
tance  [24].  Upon  careful  inspection  of  the  expression  for  the 
probe  self-impedance  Zjjf10,  it  is  noted  that  the  second  term  in 
(32)  containing  Tj™(Arp)  js  zero  when  R™  is  equal  to  one— the 
case  of  a  probe-fed  parallel  plate  waveguide.  When  R™  is  not 
equal  to  one,  this  term  diverges.  This  is  because  the  uniform 
current  on  the  probe  leads  to  a  singular  charge  accumulation  at 
the  probe  end  giving  rise  to  a  singular  reactance.  Thus,  in  the 
case  of  a  microstrip  disk  excited  by  a  probe,  in  order  to  account 
properly  for  the  probe-self  impedance,  the  continuity  of  the 
current  at  the  probe-disk  junction  must  be  ensured.  If  we  take 
the  impedance  due  to  the  probe  and  the  attachment  mode  current 
together,  we  arrive  at  the  following: 


3U  *’  2  it, 

[l  -  R]?J][l  -  /t^][l  - 
[l  - 

•  e'*lr<ri-4iy*2**2  (30} 

f)  '  ~  2  *3  (1  -  fl™***3*"3]  W 

and  ij,  =  yjuf/tt .  The  expressions  for  the  generalized  reflec¬ 
tion  coefficients  R^,  and  R £,/  are  given  in  the  Appendix.  Z^T> 
represents  the  self-impedance  of  the  probe  and  may  be  expressed 


-  J k3h3J0(k3R)H"'(k3R )  -±-f 


■dk'k'Jtik'R)^  J^*,)  (32) 

Z^4"1  is  the  input  impedance  term  due  to  the  attachment  mode 
and  is  given  by 

!~dk*k'J 
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where  the  divergent  term  in  the  probe  self-impedance  Zjj*10  has 
been  cancelled  by  the  contribution  of  the  attachment  mode  which 
ensures  continuity  of  the  current. 

V.  Radiation  Fields 

The  radiation  field,  or  far  field,  components  in  region  0  may 
be  obtained  from  the  longitudinal  components  with 

_  VqHoz  -  E0z 


For  large  observation  distances,  the  expressions  for  the  field 
components  may  be  evaluated  using  the  saddle  point  method 
with  the  saddle  point  being  kp  =  k0  sin  6  where  8  = 
tan_1(p/z).  The  longitudinal  field  components  due  to  disk 
current  Kj(p)  are  given  by 

E$(p,z)=  £  e'm4( -i)m{ kpe0J(kfi , zj) 

m  -  -  od 

eik0r 


The  first  term  of  the  probe  self-impedance  in  (32)  corresponds 
to  the  input  impedance  of  a  coaxial  probe  driven  parallel-plate 
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Likewise,  the  radiation  field  component  due  to  the  probe  is  given 
by 


#,(*.*)-  E  eim*(~iY 
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and  R  *  0.048 a2.  Given  these  parameters,  die  antenna  is  char¬ 
acterized  by  varying  the  upper  radius  a,  and  the  separation 
between  the  disks  h2-  Convergent  results  for  the  input  impedance 
and  radiation  fields  using  cavity  mode  basis  functions  are  ob¬ 
tained  with  (N  =  4 ,P  =  3),  (R  -  4 ,S  =  3),  while  those  using 
the  Chebyshev  polynomial  basis  functions  with  the  edge  condi¬ 
tion  use  (N  *  3 ,P  =  3),  (R  =  3,5  =  3).  For  the  calculation  of 
the  input  impedance  given  by  (29),  very  good  results  are  ob¬ 
tained  with  m  =  ±1  for  the  terms  associated  with  the  disk 
current  amplitudes,  b^p,  d%r,  and  and  taking  m  * 
0,  ±  1 ,  ±  2  for  the  probe  self-impedance  and  attachment  mode 
terms  in  (34).  Additional  modes  produce  only  a  1  or  2  0 
difference  in  the  input  impedance  calculations  for  the  parameters 
considered  here.  It  is  found  numerically  that  the  probe  self-im¬ 
pedance  and  attachment  mode  impedance  terms  taken  together  in 
(34)  give  rise  to  a  primarily  inductive  reactance  contribution. 
Computation  time  for  the  input  impedance  of  the  stacked  struc¬ 
ture  is  approximately  a  half-hour  of  CPU  time  per  frequency  on 
a  VAXstadon  3500. 

Calculated  and  measured  [6]  reflection  coefficients,  T  =  (Zm 
-Z0)/(Zm  +  Z0)  where  Z0  =  50  0,  are  shown  in  Figs.  2(a) 
and  2(b)  for  the  stacked  configuration  case  A  with  a,  j  a2  =  1.01 
and  h2  /a2  equal  to  0.36  and  0.48,  respectively.  The  agreement 
between  the  measured  and  calculated  results  is  very  good.  The 
loop  in  the  impedance  locus  of  Fig.  2(a)  reduces  in  size  in  Fig. 
2(b)  as  h2 /a2  is  increased  from  0.36  to  0.48,  leading  to  the 
wide  bandwidth  behavior  of  this  configuration.  Especially  off 
resonance,  it  is  seen  that  the  attachment  mode  and  probe  self-im¬ 
pedance  terms  are  required  for  accurate  results. 

In  Fig.  3,  return  loss  calculations  using  both  cavity  mode  and 
Chebyshev  polynomial  basis  functions  are  compared  with  mea¬ 
sured  results  [6]  for  the  case  of  Fig.  2(b).  The  agreement 
between  the  calculated  and  measured  results  is  very  good.  A 
frequency  shift  in  the  results  on  the  order  of  1-2%  is  observed 
for  the  two  sets  of  basis  functions.  Due  to  the  presence  of  the 
upper  disk  in  the  stacked  configurations,  two  resonances  associ¬ 
ated  with  the  two  constitutive  resonators  of  the  stacked  structure 
[11]  are  easily  distinguished,  giving  rise  to  a  16  %  —15  dB 
bandwidth  in  this  case.  One  resonance  is  associated  with  the 
resonator  formed  by  the  lower  disk  and  the  ground  plane  and  the 
second  resonance  is  associated  with  the  resonator  formed  by  the 
two  disks.  Comparing  the  return  loss  of  the  stacked  configura- 
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VI.  Numerical  Results  and  Discussion 

The  integrals  of  the  matrix  elements,  (26)  and  (27),  and  in  the 
impedance  expression  (29),  are  evaluated  numerically  along  an 
integration  path  deformed  below  the  real  axis  to  avoid  the 
singularities  on  the  real  axis  which  correspond  to  the  radiating 
and  guided  modes  of  the  layered  medium.  When  using  cavity 
mode  basis  functions,  the  integrands  vary  asymptotically  as 
1  jk\  while  those  using  the  Chebyshev  polynomial  basis  ftinc- 
tions  with  the  edge  condition  vary  asymptotically  as  1/fc2.  To 
enhance  the  convergence  of  the  integrals  when  using  the  Cheby¬ 
shev  polynomial  basis  functions,  the  asymptotic  values  of  the 
integrands  are  subtracted  out  and  evaluated  analytically. 

For  the  stacked  microstrip  configurations  discussed  here,  the 
following  parameters  are  used:  a2  *  1.3233  cm,  h3  =  2 /»,  ■ 
0.115j2,  «,  *  «j  *  2.45e0,  t2  *  1 .22 e0  (foam),  p0  *  0.6a2, 


tion  to  that  of  the  single  disk,  when  the  upper  disk  and  substrate 
are  removed,  it  is  shown  that  the  input  impedance  of  the  single 
disk  presents  an  impedance  mismatch.  When  the  probe  position 
is  changed  from  0.6a2  to  approximately  0.3a2  to  obtain  a 
match,  a  2.3  %— 10  dB  bandwidth  is  achieved. 

Fig.  4  illustrates  the  effect  of  the  separation  h2  /  a2  on  the 
input  impedance  for  the  cases  a, /e2  =  1.05,  where  the  im¬ 
pedance  has  been  calculated  using  the  Chebyshev  polynomial 
basis  functions  including  the  probe  self-impedance  and  attach¬ 
ment  mode  terms.  While  the  position  of  the  lower  resonance 
remains  essentially  the  same,  the  position  of  the  upper  resonance 
is  a  function  of  the  height  h2,  decreasing  with  increasing  h2.  As 
seen  in  the  figure,  the  excitation  of  the  upper  resonance  in¬ 
creases  with  increasing  h2  (up  to  a  certain  h2  beyond  which 
there  is  little  coupling  [11]).  This  is  due  to  the  fact  that  the 
coupling  interaction  between  the  two  modes  increases  as  the 
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0.1  GHz  Sncrtmant  3.2  -  3.0  GHz 
o,/Oj  •  1.01  h,/o,  »  0.36 

o  (3,3)  Ch*by**«v  w»>  attochmonl  modi  (m«0.±1.±2) 

a  (4,3)  Cavity  modii  (m»0.*1.±2)  with  cttoch.  mode  (m-0.*1.±2) 

•  (3.3)  Ch«by*h*v  (m-*1)  •ithout  attochmml  mode 

•  totasumd  [6] 

(a) 


0.1  CHz  iweromont  3.2  -  3.0  GHz 
0,/a,  •  1.01  hj/oj  ■  0.46 

0  (3.3)  Ch«byth«v  (m«*  1)  with  ottochmwnt  mod*  (m«0.*1.*2) 
a  (4,3)  Cavity  mod*»  (m-0.±1.*2)  with  attach,  mode  (m-0.±1.*2) 

•  (3.3)  Chabyahav  (m-*1)  without  attachment  mode 

•  Moomurvd  [6] 

<b) 

Fig.  2.  T  of  sucked  configuration  A.  a,  /a2  »  1.01.  (a)  h2  /a2  «  0.56. 
(b)  h2 /a2  »  0.48. 


upper  resonant  frequency  approaches  the  lower  resonant  fre¬ 
quency,  i.e.,  as  the  upper  resonant  frequency  decreases  with 
increasing  h2.  Or,  conversely,  the  coupling  interaction  de¬ 
creases  as  the  separation  between  the  disks  approaches  zero. 
Calculated  and  measured  [6]  return  loss  results  are  compared  in 
Fig.  5,  where  a  13%  —  10  dB  bandwidth  is  obtained  in  Fig.  5(a) 
and  where  dual  frequency  operation  is  observed  in  Fig.  5(b). 
Again,  the  agreement  between  the  calculated  and  measured 
results  is  good. 


Fr»q ueney  (GHz) 

Fig.  3.  Return  loss  of  sucked  configuration  A.  o,  /a2  -  1.01  and  h2  / a2 
=  0.48.  Return  loss  of  single  disk  with  no  upper  substrate  and  p0/a2  = 
0.3, 0.6. 


3.5  4  4.5  5 

Frequency  (GHz) 

Fig.  4.  Input  impedance  of  sucked  configuration  A.  a,  /a2  =  1.05  and 
h2  /a2  -  0.24  (-),  0.36  (-  -).  0.48  (— ). 


Shown  in  Fig.  6  are  the  input  impedance  results  of  configura¬ 
tions  A  and  B  with  parameters  a,  /a2  -  1.2  and  h2ja2-  0.24. 
Generally,  the  two  configurations  have  similar  characteristics. 
For  configuration  B,  the  increased  distance  between  the  two 
disks  and  the  higher  “effective”  dielectric  constant  between  the 
disks  results  in  an  upper  resonance  occurring  at  a  lower  fre¬ 
quency  as  compared  with  configuration  A. 

Illustrated  in  Fig.  7  are  the  disk  current  distributions  for 
configuration  A  with  a,  /a2  =  1.01  and  h2 /a2  =  0.48  at  the 
lower  resonance,  that  is  with  kia2  =  1.655  using  (4,4)  Cheby- 
shev  basis  functions  for  each  disk  and  k2a2  =  1.68  using  (5,4) 
cavity  mode  basis  functions  for  each  disk.  As  the  number  of 
cavity  mode  basis  functions  is  increased,  the  singular  behavior  at 
the  disk  edge  of  the  $  component  of  the  current  distribution  is 
better  characterized.  The  magnitude  of  the  $  component  of  the 
current  for  the  upper  disk  is  approximately  uniform  across  the 
disk  where  the  amplitude  slightly  increases  toward  the  edges  due 
to  the  parasitic  effect  of  the  upper  disk  excited  by  the  fringing 
fields. 

In  Fig.  8,  the  radiation  patterns  of  the  stacked  microstrip 
antenna  configuration  of  Fig.  7  are  compared  with  those  of  the 
single  disk  with  no  upper  substrate.  For  the  probe-fed  single 
microstrip  disk,  the  probe  position  is  taken  to  be  p0/c2  =  0.3, 
while  p0  ja2  =  0.6  for  the  stacked  configuration.  The  E4  com¬ 
ponent  remains  essentially  the  same  for  both  the  single  disk  and 
stacked  configuration.  The  radiation  pattern  of  the  stacked  con- 
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(C) 

Fig.  5.  Return  loss  of  stacked  configuration  A.  a,  /a2  »  1.05.  (a)  h2  / b2 
«  0.48.  (b)  h2  /aj  *  0.36.  (c)  h2  /a2  =  0.24. 


Fig.  6.  Inpw  impedance  of  stacked  configurations  A  and  B.  ax  /  a2  «  1.2 
and  h2/a2*  0.24. 

figuration  is  more  directive  than  (hat  of  the  single  disk,  where 
the  Ef  beamwidth  is  decreased  in  the  stacked  case. 

VII.  Conclusion 

The  input  impedance  of  a  microstrip  antenna  consisting  of  two 
circular  microstrip  disks  in  a  sucked  configuration  driven  by  a 


(a) 


Fig.  7.  Disk  currents  of  stacked  configuration  A.  a,  /fl2  =  1.01,  h2  fa2  *= 
0.48.  (a)  a,  1  KJt  // 1 .  (b)  a2  |  KJt  // 1 . 


Fig.  8.  Radiation  pattern  for  stacked  configuration  A  (— )  with  a,  /a2  « 
1.01,  h2  /a2  *  0.48,  *3a2  »  1.68,  and  for  single  disk  ( •  •  • )  with  k2a2  « 
1.75  and  p0/e2  “  0.3. 


coaxial  probe  is  investigated.  A  rigorous  analysis  is  performed 
using  a  dyadic  Green’s  function  formulation  where  the  mixed 
boundary  value  problem  is  reduced  to  a  set  of  coupled  vector 
integral  equations  using  the  vector  Hankel  transform.  Galerkin’s 
method  is  employed  in  the  spectral  domain  with  an  additional 
term  used  in  the  current  expansion  to  account  for  the  singular 
nature  of  the  current  in  the  vicinity  of  the  probe,  ensure  continu¬ 
ity  of  the  current,  and  to  speed  up  convergence  of  the  solution. 
Ensuring  continuity  of  the  current  by  means  of  the  attachment 
mode  is  shown  to  be  necessary  for  rigorously  including  the 
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pn.be  self-impedance  and  obtaining  accurate  results  for  the  nn  ^a/+u  +  ^n</+  ,)e/2*(,+  ,)z*,+1 

stacked  microstrip  configuration.  The  input  impedance  of  the  j  +  ~jp  ec»(uiu*i+i  '  1 

stacked  microstrip  antenna  is  calculated  as  a  function  of  the  ,!,+ "  n<,+  l> 

layered  parameters  and  the  ratio  of  the  two  disks.  Both  wide  where  Rff0  =  R™  =  0  and  Rff3  =  - 1  and  R™  =  1.  The 

bandwidth  and  dual  frequency  operation  are  shown.  Disk  current  Fresnel  reflection  coefficients  R%±  „  and  R]tf± „  are  defined 

distributions  and  radiation  patterns  are  also  presented.  Calcu-  by 


la  ted  results  for  the  stacked  microstrip  configuration  are  shown 
to  compare  well  with  experimental  data. 

VET  Appendix 

The  explicit  expressions  for  |™(kp.z).  £™(*P>Z,z').  and 
(JJ(kp,z,z')  of  (11)  and  (12)  are  given  here.  For  the  stacked 
configurations  of  Fig.  1 ,  we  have 

[l  ~  R™\[l  +  R™«,*3z',3][l  -  e'*31*3] 
[1  -  RLWa*'2*3*'’3] 

(40) 

[l  -  j?™<.'2*ir*i][l  ~  R™e'2*2**2] 

(41) 

where  r\,  =  y/Hi/f/  and  z’2  =  d2- 
For  source  and  observer  are  in  region  1  (assuming  z  >  z’), 
the  expressions  are 


nTE  R{i±\)kn  itika±  i>* 

~  7  7  7 

I*</±l)*/2  M/*</±  1)2 

nTM  eU±\)klz  ~  ltkU±  1)2 

(,±i)  ~  k  TT  k  • 

tU±\)Klt  *  tlKa±i)z 
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Abstract-  A  new  method  for  analyzing  frequency-dependent  transmission  line  systems 
with  nonlinear  terminations  is  presented.  The  generalized  scattering  matrix  formulation  is 
used  as  the  foundation  for  the  time  domain  iteration  scheme.  Compared  to  the  admittance 
matrix  approach  proposed  in  a  previous  paper,  it  has  the  advantage  of  shorter  impulse 
response  which  leads  to  smaller  computer  memory  requirement  and  faster  computation 
time.  Examples  of  a  microstrip  line  loaded  with  nonlinear  elements  are  given  to  illustrate 
the  efficiency  of  this  method. 

I.  INTRODUCTION 

In  recent  years,  the  effect  of  the  interconnection  lines  on  high-speed  integrated 
circuits  has  become  more  and  more  important.  As  the  speeds  of  integrated  circuits 
increase,  the  propagation  delay  as  well  as  the  dispersion  and  loss  of  interconnection 
lines  can  no  longer  be  neglected.  Traditional  lumped  element  circuit  models  must 
be  supplemented  by  dispersive  transmission  line  models  in  order  to  account  for 
these  effects.  This  has  created  the  need  for  new  numerical  procedures  that  can  be 
easily  incorporated  into  current  CAD  tools.  To  make  matters  more  complicated, 
the  interconnection  lines  are  terminated  with  not  only  linear  elements  but  also 
nonlinear  semiconductor  devices,  such  as  diodes  and  transistors. 

Several  techniques  are  now  commonly  used  to  deal  with  nonlinear  circuit  prob¬ 
lems,  for  example,  the  direct  time  domain  approaches  [1,2],  and  the  semi-frequency 
domain  approaches,  such  as  the  harmonic  balance  [3,4]  and  the  modified  har¬ 
monic  balance  techniques  [5,6].  Semi-frequency  domain  approaches  are  useful 
for  microwave  and  millimeter  wave  integrated  circuits  but  become  impractical 
for  digital  integrated  circuits  because  of  their  wide-band  nature.  On  the  other 
hand,  frequency-dependent  parameters  often  make  it  awkward  to  apply  the  direct 
time  domain  approach  to  interconnection  line  systems.  The  time-domain  finite- 
difference  method  [7]  and  the  time-domain  method  of  moments  [8]  have  been 
proposed  to  deal  directly  with  electromagnetic  scattering  from  nonlinear  loads. 
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However,  dispersion  problems  are  absent  from  the  discussions. 

Liu  and  Tesche  [1]  developed  a  combined  time-domain  frequency-domain  treat¬ 
ment  of  antennas  and  scatterers  with  nonlinear  loads.  In  their  work,  the  transfer 
function  (impulse  response)  of  the  linear  portion  of  the  investigated  system  is  first 
obtained  through  the  frequency- domain  analysis,  and  it  is  then  used  to  solve  the 
entire  nonlinear  problem  in  the  time  domain.  Subsequent  improvements  to  this 
method  have  been  suggested  by  Djordjevic,  Sarkar  and  Harrington  [9]  through 
artificially  introducing  pairs  of  quasi-matched  passive  networks,  and  by  Canig- 
gia  [10]  and  Schutt-Aine  and  Mittra  [11],  through  macromodel  and  scattering 
parameter  analysis  based  on  a  fixed  reference  impedance. 

In  this  paper  we  shall  present  an  extended  and  more  natural  method  that 
will  completely  eliminate  the  need  for  any  artificial  networks  or  fixed  reference 
impedances.  Its  close  ties  to  the  physics  of  transmitted  and  reflected  waves  on 
transmission  lines  also  help  in  achieving  the  purpose  of  reducing  computation 
time.  The  algorithm  is  explained  in  the  next  section.  In  Section  III  the  de¬ 
tails  of  applying  our  formulations  to  frequency-dependent  transmission  lines  with 
nonlinear  loads  are  illustrated  in  the  analyses  of  a  nonlinearly-loaded  dispersive 
transmission  line. 


n.  ANALYSIS  BASED  ON  WAVE  TRANSMISSION  AND  REFLECTION 


An  arbitrary  system  of  n  dispersive  transmission  lines  can  be  represented  by  the 
following  coupled  linear  ordinary  differential  equations  in  the  frequency  domain: 

dd  (1) 

--|/|=;u,[C]-|V]  +  [G].[V] 

Treating  the  n  transmission  lines  as  a  2n-port  system,  we  can  derive  from  (1) 
the  admittance  matrix  [Y] ,  which  relates  terminal  voltages  to  terminal  currents: 

2n 

h  =  E  yj*v*  (*) 

k=l 

The  time  domain  counterparts  become  convolution  relations: 


*;(<)  = 


(3) 


where  is  the  inverse  Fourier  transform  of  Y^(w).  The  terminal  voltages 
and  currents  for  any  particular  system  can  then  be  uniquely  determined  once  the 
terminal  conditions 


\m  =  [f(vm 


(4) 


are  given.  If  all  the  terminations  contain  only  linear  elements,  we  can  solve  the 
problem  in  the  frequency  domain.  Otherwise,  iteration  procedures  are  usually 
required  for  obtaining  the  solutions.  The  analyses  presented  in  [1]  and  [9]  are  based 
on  equations  (3)  and  (4).  Although  their  approaches  are  straightforward,  there 
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exist  problems  that  could  possibly  affect  the  efficiency  of  numerical  computation. 
First  of  all,  the  parameter  yjk{t)  i*  equivalent  to  the  current  measured  at  port  j 
with  a  voltage  impulse  excitation  of  unit  amplitude  at  port  k  while  ports  other 
than  k  are  short-circuited.  Owing  to  strong  reflections  at  the  terminations,  the 
duration  of  Vjk{t)  i*  usually  long  for  slightly  lossy  transmission  line  systems, 
and  even  infinite  for  a  lossless  system.  The  long  duration  puts  great  demands  on 
computer  memory  and  execution  time  as  required  to  perform  convolution  integrals 
of  (u(f)]  and  [y(t)|  during  the  iterative  solution  of  nonlinear  equations. 

In  order  to  overcome  the  disadvantage  of  using  the  parameters  {y^},  the  au¬ 
thors  in  [9]  artificially  insert  a  pair  of  complement  ary  passive  networks  between 
the  end  of  the  transmission  lines  and  the  actual  terminations.  The  purpose  of 
that  is  to  make  the  augmented  linear  network,  which  consists  of  the  transmission 
line  portion  and  the  artificial  network  directly  connected  to  the  end  of  the  trans¬ 
mission  lines,  a  quasi-matched  linear  system  so  that  the  duration  of  the  impulse 
responses  for  the  augmented  linear  system  can  be  effectively  shortened.  However, 
the  other  artificial  network  which  is  directly  attached  to  the  original  load  will 
contain  negative  resistors  and  hence  may  render  the  numerical  solution  unstable, 
especially  when  the  transmission  lines  are  lossless. 

Instead  of  dealing  with  terminal  voltage  and  current,  we  will  analyze  the  trans¬ 
mission  line  system  from  the  viewpoint  of  voltage  waves.  We  choose  the  input 
and  output  waves  at  the  terminal  ports  of  the  transmission  lines  as  the  variables 
of  the  problem  as  shown  in  Fig.  1.  The  parameters  {2?^}  and  {Cj}  are  defined 
as  follows: 

BjM  =5  [VyW  -  z0,W  •  /,Ml  (5) 

CyH  =5  [VyM  +  Zo,(-)  •  JyM]  (6) 


Figure  1.  Linear  multi-port  network  consisting  of  transmission  lines. 
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where  Zoj(u)  =  yj Ljj(u))/Cjj(u)  is  the  frequency-dependent  characteristic  impe¬ 
dance  on  line  connected  to  port  j  .  The  linear  dispersive  transmission  line  system 
is  thus  characterized  by  a  scattering  matrix  [5y] ,  i.e., 


Cl  1 

Sll  Si2  513  •••  *^l,2n  " 

r  Bx 

c2 

521  S22  S2i  •  •  •  52,2» 

b2 

Cl 

= 

53i  532  533  •••  53)2n 

•  •  •  •  • 

Bz 

.  C2n . 

•  •  •  •  • 

•  .  *  •  • 

.  *-*2n,l  S2n,2  $2n,3  ***  S2n<2n 

.  B2n . 

or 


2n 

=  £  sjtn JMu.) 


fc=  1 


(Jfc  =  1,2,.. .  ,2n) 


(8) 


It  is  easy  to  realize  from  the  above  equations  that  for  all  j  ^  k ,  5jjt(w)  is 
equal  to  2Vj(u>)  if  all  ports  are  loaded  with  their  transmission  line  characteristic 
impedances  and  only  a  voltage  source  of  amplitude  1  is  applied  at  port  k .  This 
corresponds  to  impulse  response  or  transfer  functions  in  the  time  domain.  If  the 
coupling  between  individual  transmission  lines  is  weak,  the  system  will  be  close 
to  being  perfectly  matched.  In  this  case  we  can  conclude  that  the  inverse  Fourier 
transform  of  Sj±  ,  denoted  as  ( t ) ,  will  be  of  much  shorter  duration  than  y^(<) . 
Therefore  we  can  effectively  reduce  the  memory  required  to  store  the  past  values 
of  hjk  and  the  time  to  compute  the  convolution  integrals  in 


dr/ijJfc(<-r)6*(r) 


0) 


without  inserting  any  artificial  networks. 

We  now  have  to  solve  [6(f )]  and  (c(i)]  by  incorporating  the  nonlinear  boundary 
conditions  of  the  problem.  Specifically  at  port  j ,  equations  (5)  and  (6)  lead  to 

18] 

VyM  =  ByM  +  C»  (10) 

and 

2  CjM  =  VjM  +  ZojW  ijH  (11) 

By  taking  inverse  Fourier  transform  on  both  sides  of  (10)  and  (11),  we  obtain 
their  time  domain  counterparts: 

»i(0  =  *i(<)  +  Cj(t)  (12) 

Vj(t)  =  2 Cj{t)  -  l  dr  zoj{t  -  t)  ))  (13) 

Jo 

where 

*o,'(<)  =  [Zo ,'M] 


Our  problem  will  be  solved  in  a  time-marching  fashion.  At  any  instant  t ,  the 
iteration  procedure  is  as  follows: 

(i)  Set  up  initial  guess  of  [b(t)\ .  A  reasonable  choice  is  to  take  values  from  the 
previous  time  step. 
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(ii)  Compute  individual  {cj(<)}  using  (9). 

(iii)  Apply  standard  nonlinear  equation  techniques  such  as  Newton- Raphson  me¬ 
thod  to  (13)  to  solve  for  individual  {t>y(i)}  • 

(iv)  Obtain  the  next  guess  of  [6(f)]  from  the  relation  bj(t)  =  Vj(t )  —  cy(<) ,  and 
compare  the  guess  with  the  previous  one.  If  the  error  is  above  a  pre-set 
tolerance,  repeat  steps  (ii)  to  (iv)  with  the  new  guess. 

For  all  practical  purposes,  we  further  divide  these  variables  into  two  sets.  The 
first  corresponds  ic  the  source  side  of  the  transmission  line  system  (Fig.  2),  iden¬ 
tified  by  odd-numbered  subscripts,  and  the  other  corresponds  to  the  load  side 
(sourceless)  with  even-numbered  subscripts  (Fig.  3).  The  reason  for  doing  so  is 
based  the  following  observations.  The  transfer  functions  linking  {&2j}  to  {c2j+i} 
and  {f>2j+l }  to  (c2j)  include  time  delay  operators  in  order  to  account  for  the 
finite  speed  of  propagation.  In  other  words,  {c2j+\}  depend  not  on  .the  present 
but  the  past  values  of  {&2j}  *ad  the  like  for  {c2j}  on  {&2j+l}*  As  we  carry 
out  the  iteration  procedure  step  by  step,  the  present  values  in  one  set  will  not 
interfere  with  those  in  the  other.  Therefore,  the  update  of  variables  can  be  done 
simultaneously  for  any  given  time  if  parallel  processing  facilities  are  available. 

In  the  next  section  we  shall  present  the  application  of  this  method  to  a  single 
dispersive  line  loaded  by  nonlinear  impedance. 


m.  N ONLINE ARLY-LOADED  DISPERSIVE  TRANSMISSION  LINE 

Shown  in  Fig.  4  is  a  uniform  dispersive  transmission  line  of  length  l  driven  by 
a  source  eo(t)  with  a  linear  source  resistance  Rt  at  one  end,  and  terminated 
by  a  nonlinear  resistor  Rn  at  the  other  end.  The  transmission  line  portion  of 
this  problem  can  be  described  in  terms  of  the  frequency  dependent  characteristic 
impedance  Zq(w)  and  effective  propagation  constant  /?(u>) .  It  can  be  shown  that 
the  frequency- domain  icattcring  matrix  is  given  as  follows 
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Figure  3.  Even-numbered  ports  (load  side). 


Figure  4.  Nonlinearly  loaded  dispersive  transmission  line. 

■5n(w)  Si2(w) 

.Si  i(w)  S22(u>) 

Noting  that  both  Su  and  S22  are  zero  for  all  time;  therefore,  we  only  need 
to  consider  S\i  and  S2 1 ,  which  involve  time  delays.  The  simplicity  of  (14) 


0 

e-JPW  0 


(14) 
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can  be  attributed  to  our  definition  of  sc.  ttering  matrix.  As  an  example,  we 
consider  a  microstrip  line  of  width  W .  The  substrate  is  of  thickness  h  and 
has  a  dielectric  constant  tr .  Numerous  empirical  formulas  are  available  from 
the  literature  (see  for  example  [13-15]).  In  this  paper,  we  will  use  the  following 
expressions  to  calculate  the  frequency  dependence  of  microstrip  line  characteristic 
impedance  and  effective  dielectric  constant  «*(<*/) ,  which  is  defined  from  /?(u>)  = 


**(")  =  /  \  2 
«r  \<^ij 

(15) 

Z0(u>)  -  _ 

We(u>)y/ee(ut) 

(16) 

and  the  effective  width  We(v)  is  governed  by  the  equation 

W'M  =  w+”i°)~w 

1  +  U)  (U)g 

(H) 

where 

7rZ0(0) 

MO* 

(18) 

u9  = 


ire 


and  the  low-frequency  limits  of  ee  and  Zq  are  respectively 


e«(0)  = 


er  +  1 


+ 


- F(W/h ) 


with 


and 


F(W/h)  -  /<1  +  12VWTI/2  +0.04(1  -  W/h)1 


if  W/h  <  1 
if  W/h  >  1 


(19) 


(20) 


(21) 


Zo( 0)  = 


+  0.25 


W 


)• 


+  1.393  +  0.667  In  +  1.444 


)} 


-1 


if  W/h  <  1 


if  W/h  >  1 


(22) 

Note  that  We(0)  in  (17)  is  derived  from  (16)  by  substituting  in  Zq(0)  and  ee(0)  . 

In  this  paper,  we  assume  that  the  width  W  is  equal  to  0.508  mm  and  the  depth 
and  the  dielectric  constant  of  the  substrate  are  equal  to  0.216  mm  and  10.2  respec¬ 
tively.  It  can  be  readily  derived  from  the  formulas  that  the  effective  permittivity 
in  the  low-frequency  limit  is  equal  to  7.46,  and  the  high-frequency  characteristic 
impedance  is  Zq(oo)  =  50.2  Cl .  The  calculation  of  the  impulse  response  function 
*■12  (^)  (—  ^2l)i  the  inverse  Fourier  transform  of  Si2(u>),  however,  is  not  trivial. 
Because  lim  Si2(u>)  does  not  approach  0 ,  the  inverse  Fourier  transform  is  sin- 

Uf — *oo 

gular.  If  numerical  computations  are  not  carried  out  with  great  care,  the  accuracy 


190 


G%  et  al. 


would  be  questionable.  In  [16]  a  similar  Fourier  integral  was  calculated.  The  only 
difference  is  that  the  input  is  not  an  impulse  function,  and  hence  it  is  numerically 
integrable.  The  authors  proposed  using  Taylor  expansion  method  and  the  method 
of  stationary  phase  for  narrow-band  input  signals.  But  for  other  input  waveforms 
a  brute  force  numerical  integration  method  was  applied  which  proved  to  be  very 
time-consuming.  To  overcome  the  difficulty  of  numerical  integration  once  and  for 
all,  our  approach  is  to  separate  the  transfer  function  5 j2  into  two  parts,  one  that 
can  be  analytically  inverted  to  an  impulse  function,  and  one  that  is  well-behaved 
and  integrable.  The  latter  requires  that  the  integrand  approaches  0  as  u>  goes 
to  infinity.  This  leads  to  a  natural  way  of  separation: 

e-j0(u>)t  _  (e-j0(u,)l  _  +  e-j0ool 

=  t~^1  {  [  e-3{0(*)-0~)l  _  1  ]  +  1  }  (23) 

where  0 00  =  ^y/M£Oer  • 

In  (23),  the  factor  e-^°°*  corresponds  to  a  time  delay  of  rT  =  l/{cy/ «>)  and 
the  inverse  Fourier  transform  of  1  is  an  impulse  function.  Therefore,  defining 
h(t)  —  —  1}  ,  we  then  have 


h12(f)  =  -  rr)  +  h(t  -  Tr)  (24) 

Now  only  h(t)  needs  to  be  evaluated  numerically.  Because  the  kernel  of  this 
Fourier  inversion  has  essentially  a  finite  range  of  integration,  the  calculation  be¬ 
comes  easier.  The  impulse  response  function  for  a  10  mm  long  microstrip  line  is 
shown  in  Fig.  5.  Only  a  limited  portion  of  hi2(<)  surrounds  the  impulse  function. 
This  is  consistent  with  our  claim  that  the  impulse  response  has  a  very  limited 
duration. 

Because  the  scattering  parameters  hjj  and  h2 2  8X6  zero,  the  input-output 
wave  variable  pair  (fej(t),  Cj(t)}  are  only  linked  through  (12).  The  iterative  solu¬ 
tion  to  a  single  transmission  line  problem  is  therefore  relatively  simple.  Once  we 
finish  calculating  Vj(i) ,  iteration  step  (iv)  will  readily  return  the  correct  values 
for  bj(t)  .  There  is  no  need  to  go  back  to  step  (i). 

We  first  examine  the  response  from  a  Gaussian  pulse  with  an  amplitude  of  1.0 
volts  and  a  width  of  10  ps  measured  at  its  half  amplitude  is  used  as  the  source. 
The  source  resistance  R,  is  50  fl .  The  load  characteristics  is  fully  described  by 
the  following  equation: 

<2  = /o  [«P  (5^5)  -  l]  -A  (25) 


(volt)  htt  (volt/psec) 


Nonlinear  Load 
Linear  Load 


(volt)  Vt  (volt) 
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Length  Effect 


Figure  9.  Load  voltages  for  different  line  lengths. 


Pulse  Width  Effect 


Figure  10.  Load  voltages  for  different  pulse  widths. 
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Figure  11.  Load  voltages  for  another  type  of  nonlinear  termination. 

For  the  co.se  Iq  =  1.0,  the  time  responses  at  both  ends  are  plotted  in  Fig.  6. 
The  puls<-  has  been  broadened  and  the  negative  trailing  edge  »  rather  significant. 
In  order  to  analyze  what  comes  from  dispersion  of  the  transmission  line  and  what 
comes  from  the  load  nonlinearity,  we  compared  the  plots  of  the  load  responses 
when  7o  is  varied  from  0.2  to  50.0  in  (25),  and  when  the  load  is  a  50  ft  linear 
resistor.  The  results  are  depicted  in  Figs.  7  and  8.  Apparently  the  nonlinearity 
contributes  most  to  the  broadening  effect.-  The  negative  trailing  edge  originates 
from  dispersion,  but  is  enhanced  by  the  nonlinear  load.  Since  the  load  charac¬ 
teristics  is  similar  to  that  of  a  typical  diode,  it  behaves  like  an  open  circuit  with 
respect  to  an  incoming  negative  voltage  wave.  This  gives  rise  to  roughly  twice  the 
response  compared  to  the  one  at  the  50  ft  load,  which  is  nearly  matched  to  the 
transmission  line.  The  large  negative  trailing  edge  also  accounts  for  the  extended 
ringing  after  it  is  reflected  back  to  the  source  end. 

In  Fig.  9,  the  load  end  voltages  for  5  mm  and  20  mm  lines  are  compared  against 
the  10  mm  line  case.  The  number  of  aero  crossings  increases  with  length,  as 
expected  from  the  transmission  line  dispersion.  Similar  phenomenon  is  observed 
when  shorter  pulses  are  injected,  as  in  Fig.  10.  The  centers  of  these  pulses  are 
intentionally  separated  to  allow  clearer  comparison.  It  is  interesting  to  note  the 
similarity  of  the  5  ps  pulse  output  and  the  20  mm  line  output  in  Fig.  9. 

We  have  also  tested  our  iteration  algorithm  on  another  class  of  nonlinear  loads 
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with  t  -  v  characteristics  described  as 

i,  =  70[taoh(^^)+1.0_ 


(26) 


Unlike  the  one  described  in  (25),  there  is  cap  on  the  positive  current.  The  result  is 
that  positive  voltage  has  a  larger  amplitude,  which  is  controlled  by  Iq  ,  as  shown 
in  Fig.  11. 

The  time  resolution  used  in  all  but  the  5  psec  pulse  input  case  is  0.4  psec  with 
a  total  of  1000  points.  For  the  latter,  the  resolution  is  0.2  psec  and  there  are  2000 
points.  Because  very  limited  portion  of  the  impulse  response  hj2(0  is  significant, 
the  actual  number  of  points  involved  in  the  convolution  integral  is  far  lower.  We 
used  the  Newton- Raphson  iteration  procedure  for  the  nonlinear  equation  (13). 
On  a  VAX  Station  3500,  the  testing  cases  take  about  4  to  27  seconds  of  CPU  time 
to  generate  the  solutions  with  1000  points.  The  large  variation  in  computation 
times  is  due  to  different  convergent  rates  for  different  types  of  nonlinearity. 


IV.  CONCLUSIONS 

A  new  method  for  the  transient  analysis  of  a  frequency-dependent  transmission 
line  system  terminated  with  nonlinear  loads  has  been  presented.  This  method  is 
not  only  effective  for  saving  the  CPU  time  required  for  solving  nonlinear  transient 
problems,  but  is  also  compact  mid  natural  in  form.  Our  generalized  scattering 
matrix  approach  is  closely  tied  to  the  concept  of  waves.  Therefore,  no  extended 
reflection  arises  as  a  result  of  artificial  boundary  conditions  as  can  occur  with  the 
admittance  matrix  method,  and  duration  of  the  impulse  responses  for  the  waves 
in  the  transmission  line  system  is  very  limited.  This  is  the  key  to  reducing  the 
amount  of  computation  time  and  memory. 

The  detailed  procedure  for  solving  this  kind  of  nonlinear  transient  problem 
is  given  through  an  analysis  of  a  nonlinearly-loaded  microstrip  transmission  line 
with  linear  source  resistance.  Extension  of  this  approach  to  multiple  transmission 
line  systems  with  nonlinear  source  and  terminations  is  being  studied  and  will  be 
reported  in  future  work. 
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